Spin wave theory for antiferromagnetic XXZ spin model on a triangle lattice in the 

presence of an external magnetic field 
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Spin wave theory is applied to a quantum antiferromagnetic XXZ model on a triangle lattice in 
the presence of an in-plane magnetic field. The effect of the field is found to enhance the quantum 
fluctuation and to reduce the sublattice magnetization at the intermediate field strength in the 
anisotropic case. The possible implication to the field driven quantum phase transition from a spin 
solid to a spin liquid is discussed. 



1. Introduction 

Quantum spin systems in low dimensions continue to 
be an interesting subject. Quantum fluctuation gener- 
ally destroys the long range order in one dimension jl], 
and plays an important role in two dimension [Q] . In this 
paper, we study the ground state of spin-S" antiferromag- 
netic quantum XXZ model on a triangle lattice in the 
presence of an external in-plane magnetic field. 

Our main focus is on the effect of the magnetic field 
to the quantum fluctuation. It is well known that a spin 
system becomes polarized at a very high field (above a 
threshold), which suppresses the quantum fluctuation. 
It is, however, less clear if the magnetic field enhances 
or suppresses the quantum fluctuation at intermediate 
field strength. The latter problem may be important 
to systems near the boundary of a quantum solid (or- 
dered phase) to quantum liquid (disorder phase) tran- 
sition. Our interest on this problem is partially moti- 
vated by the recent experiments on Cs^CuCli. That 
system || is a quasi 2-dimensional 5=1/2 frustrated 
Heisenberg antiferromagnet with a weak anisotropy fa- 
voring spins aligned in the basal plane, likely due to the 
Dzyaloshinskii- Moriya interaction . Coldea et al. || |(| 
have used neutron scattering to study the ground state 
and the dynamics of the system in high magnetic field. 
Among the observations, these authors found that the 
system undergoes phase transitions from a spin solid to a 
spin liquid to a spin fully polarized states as the in-plane 
magnetic field increases. Their experiment has raised an 
interesting theoretical question of the effect of the field on 
the quantum fluctuation and on the possible field driven 
quantum spin liquid. The basal plane of Cs^CuCl^ con- 
sists of a triangle lattice with the spin-spin coupling much 
stronger along one direction than along other two direc- 
tions. There have been some theoretical works attempt- 
ing to address some of the issues relevant to the experi- 
ments @, ^ ■ Here we shall consider a simplified model 
to study the field effect on the quantum spin fluctuation. 
For simplicity, in our model we shall neglect the lattice 
anisotropy in the real systems, and assume the spin cou- 
plings are the same between all the pairs of the nearest 
neighboring sites. We shall use XXZ model to describe 



the anisotropy in spin space, which distinguishes the in- 
plane field from the perpendicular field. 

We apply spin-wave theory |n| to study the leading 
order correction to the sublattice magnetization due to 
the quantum fluctuation. The spin wave theory or equiv- 
alent theories have been applied to a number of systems, 
including the antiferromagnetic Heisenberg model in a 
triangle lattice JllJ and the quantum XY model 1 12 . We 



calculate the ground state energy and the sublattice mag- 
netization as functions of the magnetic field, and find 
that the quantum fluctuation measured by the reduction 
of the sublattice magnetization is enhanced as the field 
increases at the intermediate field strength and in the 
anisotropic case. The possible implication to the field 
driven quantum phase transition will be discussed. 

This paper is organized as follows. In section 2, we 
introduce the model Hamiltonian. Section 3 describes 
the application of the spin wave theory to the present 
model. The results for the ground state energy and the 
sublattice magnetization are presented in section 4. In 
the final section, we summarize our results and discuss 
their possible implications to the field driven quantum 
spin liquid. 



2. Model Hamiltonian 

We consider spin-5 antiferromagnetic XXZ model on 
a triangular lattice of 3iV sites in the x-z plane in the 
presence of an external magnetic field h along the z-axis. 
The lattice constant is set to be 1. The Hamiltonian can 
be written as follows: 



H — ^~^(Sf Sf + g + SfS;, 



i+S 



(1) 



L6 



In the above equation, is the /j, component of the 
spin-S operator on site i, and S denotes the three nearest- 
neighbor bond vectors (see Fig.l) given by, S a = (1,0), 
Sp = (-1/2, V3/2), and <5 7 = (-1/2,-73/2). hS is the 
strength of the magnetic field, and h will be considered 
to be order of unity for the purpose of analyzing the spin 
wave theory, so that the magnetic field contribution to 
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FIG. 1: Three sublattices in a triangle lattice and the classical 
spin state studied in this paper. 



the energy is of the same order as the leading order spin- 
spin interaction energy, w is a parameter describing the 
anisotropy of the model, and < w < 1. The model 
is reduced to the Heisenberg model at w = 1, and to 
the XY model at w = 0. Note that the mdoel with 
w > 1 belongs to a different universal class and will not 
be discussed here. 



The angle 9 is determined variationally by dE c i/d9 = 0, 
from which we find 



cos6> = 



(>-3)/6, h < h c = 9 



1, h > h c 



(3) 



Note that the value of 9 is independent of the spin 
anisotropic parameter w for the y-component of the spin 
does not play any role in the classical limit. At zero 
field, we obtain 8 — 120°, recovering the well known re- 
sult. The case of 9 = (h > h c ) corresponds to the 
ferromagnetic phase where all the spins are fully polar- 
ized along the field direction. We see that as the field 
increases, 9 decreases from 120° to 0, and the classical 
antiferromagnetic ground state is gradually transformed 
into a ferromagnetic state. 

The leading order quantum fluctuation to the classical 
solution is proportional to S, and it is given by, 



H FL = SSY.DlM^Dj: + E a FL (h), 



(4) 



where 



-9NS/2, h < h c 
E° FL (h) = { (5) 
-3NS(h/2 - 3), h> h c 



3. Spin-wave theory 

We apply the spin wave theory to study this model [jl0| . 
This requires a proper choice of the classical state upon 
which the spin wave represents minimum quantum fluc- 
tuations. Since w < 1, we choose quantization axes 
within the x — z plane, and consider the spin quantiza- 
tions in three sublattices as indicated in Fig. 1. Namely, 
the spins in sublattice A are aligned along the z-axis and 
the spins in sublattices B and C are aligned along the 
9 and —8 directions respectively, where 9 is determined 
variationally. Such a quantization is consistent with the 
expected three sublattice classical ground state. 

We introduce three boson annihilation operators a, b, c 
and their corresponding creation operators by means of 
the Holstein-Primakoff transformation 1 1 3 to represent 
spins on sublattices A, B, and C, respectively. They are 
given by Sf = S — rii, = ^/2S — rii di, and — 
d\\/2S — rii, where d — a,b,c for the site i 6 A,B,C 
respectively, rij = d\di, and z 1 is the corresponding spin 
quantization axis, and = Sf ± iSf . We apply the 
standard spin wave theory and write the Hamiltonian in 
terms of large S'-expansion. The leading order term is 
proportional to S 2 , which is the energy of the classical 
ground state given by 

E cl = TVS 2 [6 cos 6» + 3 cos 26» + 2 cos A)]. (2) 



In Eq. (4), the sum is over the reduced Brillouin zone, 
and Dl = (At,A_£), At = (at,6t,ct), and A is the 
transpose of A. M. is a 6 x 6 matrix, and can be written 
in terms of an idendity matrix X and the Pauli matrix a x 
as 



M = T® (Mq + M + ) + . 



M_ 



(G) 



where Mo is a 3 x 3 diagonal matrix, whose matrix el- 
ements are given by (M )n = f — cos9, and {M ) 2 2 = 
(Mo)33 = ^(hcos9 — 3cos9 — 3COS26 1 ), and 





M± = | (cos6»±u;)7i 

(cos 9 ± w)7£ (cos 29 ± w)^i 



(cos 9 ± w)-f£ (cos9±w)jt 



(cos29±w)^ j(7) 




In the above equations, 7^ = (1/12) exp (ifc • <5), 
where the sum of 5 runs over {8 a , S/3, S-y}- The off diag- 
onal matrices M± represents the interaction among the 
three sublattices. The 9 term in the matrix elements 
arises from the x- and z- components of the spin fluctua- 
tion and the w term arises from the y- component of the 
spin fluctuation. 

Hfl can be diagonalized by using the Bogoliubov 
transformation. We introduce three bosonic operators 
a r with n = 1,2,3, and vector operators AL — 



lated to the original Holstein-Primikoff boson operators, 
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Dg, by a Bogoliubov transformation U, and that Hfl is 
diagonalized in terms of these new boson operators, 

_ D~ k = UV n (8) 

H FL = iS^n^l^n* + 1/2) + E° FL (h) (9) 

where w n j (n = 1, 2, 3) are the energy dispersions of the 
three spin wave excitations (magnon modes) in the spin 
ordered ground state. 

In general, the Bogoliubov transformation can be car- 
ried out hence the magnon dispersions can be calculated 
numerically. This will be discussed in the following sec- 
tion. In the ferromagnetic phase, h > h c , we have ana- 
lytical results, 

<"n,S = \/H5+!+ 3(1 + ^TJ 2 [3(1 - w)T n tf 

where T n g = 7g£™ + 7g£ 2n , with £ = exp (i27r/3). Note 
that —1/4 < T n j: < 1/2. This result is in agreement 
with the previous work Jul. 
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FIG. 2: (a) The classical energy per bond as function of mag- 
netic field h (b) The leading order quantum correction to the 
energy in spin wave theory as functions of h for several values 
of w. 



4. Results 

In this section, we present the results for the ground 
state energy and the sublattice magnetization of the an- 
tiferromagnetic XXZ model obtained from the spin wave 
theory described in the previous section. 



4.1. Ground state energy 

In this subsection, we discuss the ground state energy 
up to the leading order quantum fluctuation, namely to 
the order of S in the large S'-expansions. The ground 
state energy per bond (there are total 9N bonds in the 
lattice) is given by 

£o = e c ; + tfi (10) 

where e c i and e/j are the classical and the fluctuation 
energies, respectively. From Eq. (2), e c ; is given by, 

cl l(l-f)S 2 , h>h c [ ' 

To obtain en, we include the zero point energy of the 
bosons in Eq. (9), which leads to 

r (4+9sEv)s. h < h c 
€fl = (-(l-i)+ fe ^E^. h>h c ^ 

\ n,k 



We have solved the Bogoliubov transfromation problem 
numerically to diagonalize Hf l and to calculate the dis- 
persions u> n g and the energy eo. In Fig. 2(a), we plot the 
energy of the classical ground state as a function of the 
magnetic field h. The energy decreases monotonically as 
h increases, and is linear in h in the ferromagnetic state 
with h > he- In Fig. 2(b), we show the quantum cor- 
rection to the ground state energy as functions of h for 
several values of w. Firstly, we note that the quantum 
correction to the energy is very small in comparison with 
the classical energy even for S — 1/2. This indicates 
that the spin wave theory is a good approximation for 
the present model. At h = 0, the quantum correction is 
largest in the Heisenberg model (w = 1) and smallest in 
the XY model (w = 0). This is consistent with the intu- 
ition that the y— component of the spin would increase 
the quantum fluctuation. As h increases, however, the 
magnitude of the quantum fluctuation as a function of 
w reverses the order as we can see from Fig. 2(b). This 
becomes apparent at h = he, where the quantum cor- 
rection to the energy in the Heisenberg model vanishes, 
while that in the XY model reaches the maximum. It 
is interesting to note that the magnetic field enhances 
the quantum fluctuation in certain parameter region of 
the mdoel, a point we will come back for further discus- 
sion in the next subsection. It is also interesting to note 
that the leading order quantum fluctuation remains fi- 
nite in the general case even in the ferromagnetic phase. 
The Heisenberg limit of the model is only a special case 
where the leading order quantum fluctuation vanishes. 
The latter can be understood by examining the matri- 
ces in Hfl- At h > he, and w = 1, we have M_ = 0, 
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so that boson annihilation operators (a, b, c) do not mix 
with their creation operators (a',w,c'), hence the sys- 
tem is similar to the ferromagnetic Heisenberg model and 
has zero quantum fluctuation. 



4.2. Sublattice magnetization 

We now move to our main interest to discuss the quan- 
tum corrections to the sublattice magnetization. The 
magnetization in sublattice L — {A, B, C} is defined as 
the average spin component within the same sublattice 
along its quantization axis, 



(Sl) = jj52(S?) = S-(*Si 



(13) 



where S is the classical value of the spin, and the second 
term, ASl = (1/JV) (<4<i)> is the reduction from the 

classical spin value due to the quantum fluctuation, and 
(Q) is the expectation value of operator Q in the ground 
state. ASl can be calculated from the Bogoliubov trans- 
formation matrix U, 



ASr 



\Ul 5 



\Uw\ 2 



(14) 



where / = {1,2,3} for L = {A, B,C} respectively, and 
U is obtained in the diagonalization of Hfl ■ The results 
are shown in Fig. 3 for ASa and in Fig. 4 for ASb 
for functions of the field h for several values of w. By 
symmetry, ASc — ASb- 

We first discuss the Heisenberg limit of w = 1. At 
h = 0, we find ASl ~ 0.26 for all the three sublattices. 
This result is the same as that reported early using a 
Schwinger boson mean field theory by Yoshioka As 
h increases, both ASa and ASb decrease, but ASa drops 
much faster as we can see in the figure. For h > he = 9, 
the ground state in ferromagnetic, and we find ASa = 
ASb = as a result of the zero quantum fluctuation in 
this case as we discussed in the previous subsection. 

We now turn to discuss the XY limit of the model with 
w = 0. At h = 0, AS A = AS B « 0.05, which is much 
smaller than the value in the Heisenberg model (~ 0.26), 
but is comparable to the value of 0.06 reported in the 
square lattice XY model |12| ]. As h increases from 0, 
ASa decreases to zero while ASb increases to reach a 
maximum at h = 3. Figs 3 and 4 also show that ASl is 
symmetric with respect to the value of h = 3 in the region 
< h < 6. In the XY limit, the quantum fluctuation 
arising from the y-component of spin is absent. The off- 
diagonal matrices M± in Hfl are given by for h < he, 



I 



M± 



6 'k 



h-3 * 
6 >k 

6 <fc 



h-3 „,» 
6„ 'jfe 

r (ft-3) 

i 18 



lb* (15) 




FIG. 3: Quantum correction to local spin in sublattice A as 
function of magnetic field h for several w. 



At h = 3, there is no finite matrix element between the 
sublattice A and sublattices B or C. Therefore, the trian- 
gle lattice is decomposed into a honeycomb lattice con- 
sisting of sublattices B and C and N isolated spins of 
sublattice A. Consequently, the spins in A are all paral- 
lel to the field with ASa — 0, and ASb is the same as the 
result in a honeycomb lattice in the absence of external 
fields. The form of M± also indicate a symmetry with 
resepct to h = 3, although the rigorous mathematics in- 
volves more delicated analyses since the matrices contain 
elements with both even and odd functions of h — 3. 

As shown clearly in the figures, there are sharp peaks 
for both ASa and ASb at h = h c except for w = I. The 
peak height increases as w decreases, and is the largest for 
w — 0. Part of this feature found in our calculations may 
be understood by examining the ferromagnetic phase for 
h > he- In this phase, we have ASa = ASb = AS, 
which is given by 



AS* 



I 

6N 



E 



-3 + | 



■ 3(1 + w)T % 



(16) 



n,k 



r (h-3f 
1 18 



From this expression, we see that AS decreases as 
h further increases as we can see from the figures. 
{dAS/dh)\ h ^ hc+0 + — * — oo, so that the slopes of AS 
divergent. This explains the singularities of AS found in 
our numerical calculation shown in the figures. 

The most interesting result we can learn from these cal- 
culations is the possible field enhanced quantum fluctua- 
tions in certain quantum spin systems. This is somewhat 
in contrary to the general intuition since the field is ex- 
pected to partially polarize the spins hence to suppress 
the quantum fluctuation. Our calculations and analy- 
ses demonstrate that this intuitive argument may break 
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FIG. 4: Same as in Fig. 3 in sublattice B. 



down in certain parameter region. This finding should 
have important implications for the magnetic field driven 
quantum phase transition from a spin solid to a spin liq- 
uid. In the spin wave theory, AS* is a measure of the 
quantum fluctuation. If AS — S, the sublatitce magne- 
tization vanishes, indicating the melt-down of a quantum 
solid. 



interesting, and it may have important implications for 
systems near the solid-liquid phase boundary. Near that 
boundary, the external magnetic field may well drive a 
quantum solid to a quantum liquid due to the increase of 
the quantum fluctuation. In quasi 2-dimensional triangle 
lattice of CsiCuCl^ the in-plane spin-spin couplings are 
strongly anisotropic, and the coupling along one direc- 
tion is about three times as large as those along the other 
two directions. The spin system is between one and two 
dimensions, and the system in the absence of field has a 
long range order but is closer to the spin liquid boundary. 
In that type of systems, an enhanced quantum fluctua- 
tion may easily destroy the long range order. We spec- 
ulate that the in-plane field induced phase transition |j] 
to the quantum liquid state observed in CsiCuCl^ could 
be due to the increase of the quantum fluctuation. More 
detailed calculations on more realistic model would be 
interesting to carry out to see if such a conjecture is rel- 
evant to the experiment. 
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discussions on their experiments. We also wish to thank 
S. Q. Shen for many useful discussions on this subject, 
especially in the early stage of the work. This work is in 
part supported by Chinese Academy of Sicences, and by 
the US DOE grant DE/FG03-01EG45687. 



5. Discussions and Summary 

We have applied the spin wave theory to study quan- 
tum XXZ model on a triangle lattice in the presence 
of an in-plane magnetic field. Our model includes the 
Heisenberg model at one limit and the XY model at the 
other limit. We have calculated the ground state en- 
ergy and the sublattice magnetization as functions of the 
external field and the anisotropy of the spin coupling. 
We have found that the field may enhance or suppress 
the quantum fluctuation and the reduction of the sub- 
lattice magnetization, depending on the spin anisotropy 
and the field strength. At the intermediate field and in 
the anisotropic case, the field enhances the quantum fluc- 
tuation. The reduction of the sublattice magnetization 
in the model we have studied is still too small to destroy 
the spin long range order even for the smallest spin sys- 
tem with S = 1/2. This is because the model we studied 
here is still quite far from the solid-liquid boundary in the 
absence of the field. Nevertheless, the qualitative effect 
of the field on the quantum fluctuation we have found is 
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